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Cold dark matter (CDM) models predict small-scale structure in excess of observations of the cores 
and abundance of dwarf galaxies. These problems might be solved, and the virtues of CDM models 
retained, even without postulating ad hoc dark matter particle or field interactions, if the dark 
matter is composed of ultra-light scalar particles (m ~ 10~ 22 eV), initially in a (cold) Bose-Einstein 
condensate, similar to axion dark matter models. The wave properties of the dark matter stabilize 
gravitational collapse providing halo cores and sharply suppressing small-scale linear power. 
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Introduction. — Recently, the small-scale shortcomings of 
the otherwise widely successful cold dark matter (CDM) 
models for structure formation have received much at- 
tention (see |l|-H] and references therein). CDM models 
predict cuspy dark matter halo profiles and an abundance 
of low mass halos not seen in the rotation curves and lo- 
cal population of dwarf galaxies respectively. Though the 
significance of the discrepancies is still disputed and so- 
lutions involving astrophysical processes in the baryonic 
gas may still be possible (e.g. ||), recent attention has 
mostly focused on solutions involving the dark matter 
sector. 

In the simplest modification, warm dark matter (m ~ 
keV) replaces CDM and suppresses small-scale struc- 
ture by free-streaming out of potential wells || , but this 
modification may adversely affect structure at somewhat 
larger scales. Small-scale power could be suppressed 
more cleanly in the initial fluctuations, perhaps originat- 
ing from a kink in the inflaton potential Q , but its regen- 
eration through non-linear gravitational collapse would 
likely still produce halo cusps M. 

More radical suggestions include strong self- 
interactions either between dark matter particles ffl or 
in the potential of axion-like scalar field dark matter j| . 
While interesting, these solutions require self-interactions 
wildly in excess of those expected for weakly interacting 
massive particles or axions respectively. 

In this Letter, we propose a solution involving free par- 
ticles only. The catch is that the particles must be ex- 
traordinarily light (m ~ 10~ 22 eV) so that their wave 
nature is manifest on astrophysical scales. Under this 
proposal, dark matter halos are stable on small scales for 
the same reason that the hydrogen atom is stable: the 
uncertainty principle in wave mechanics. We call this 
dark matter candidate fuzzy cold dark matter (FCDM). 

Equations of Motion. — It is well known that if the dark 
matter is composed of ultra-light scalar particles m -C 
leV, the occupation numbers in galactic halos are so high 
that the dark matter behaves as a classical field obeying 
the wave equation 



US = mf 
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the Compton wavelength mT x but much smaller than 
the particle horizon, one can employ a Newtonian ap- 
proximation to the gravitational interaction embedded in 
the covariant derivatives of the field equation and a non- 
relativistic approximation to the dispersion relation. It is 
then convenient to define the wavefunction ip = Ae %a , out 
of the amplitude and phase of the field 6 = Acos(mt—a), 
which obeys 
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where \l/ is the Newtonian gravitational potential. For 
the unperturbed background, the right hand side van- 
ishes and the energy density in the field, p = m 2 \ip\ 2 /2, 
redshifts like matter p cx a~ 3 . 

On time scales short compared with the expansion 
time, the evolution equations become 



1 



id t il> = { V 2 + m^U, V 2 ^ = 4ttGS P . 
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Assuming the dark matter also dominates the energy 
density, we have Sp = m 2 5\ip\ 2 /2. This is simply the non- 
linear Schrodinger equation for a self-gravitating particle 
in a potential well. In the particle description, ij) is pro- 
portional to the wavefunction of each particle in the con- 
densate. 

Jeans / de Broglie Scale. — The usual Jeans analysis tells 
us that when gravity dominates there exists a growing 
mode e 7 ' where 7 2 = 4nGp; however a free field oscillates 
as e~ lEt or 7 2 = -(fc 2 /2m) 2 . In fact, 7 2 = AirGp - 
(k 2 /2m) 2 and therefore there is a Jeans scale 



where we have set h = c = 1. On scales much larger than 



rj = 2n/kj = n^HGp)- 1/4 m-^ 2 , 

= 55m2 2 1/2 (pM)- 1/4 (n m /i 2 )- 1/4 kpc, (4) 

below which perturbations are stable and above which 
they behave as ordinary CDM. Here 77122 = m/10 _22 eV 
and pb — 2.8 x 10 n r2 m ft 2 Mo Mpc -3 is the background 
density. The Jeans scale is the geometric mean between 
the dynamical scale and the Compton scale (c.f. 0-3|) as 
originally shown in a more convoluted manner by [1C|. 

The existence of the Jeans scale has a natural interpre- 
tation: it is the de Broglie wavelength of the ground state 



1 



of a particle in the potential well. To see this, note that 
the velocity scales as v ~ (Gp) 1 ^ 2 r so that the de Broglie 
wavelength A ~ (m») _1 ~ m" - 1 (Gp)- 1 /' 2 r- 1 . Setting 
r j = A = r, returns the Jeans scale. Stability below the 
Jeans scale is thus guaranteed by the uncertainty prin- 
ciple: an increase in momentum opposes any attempt to 
confine the particle further. 

The physical scale depends weakly on the density, but 
in a dark matter halo p will be much larger than the 
background density p\>. Consider the density profile of a 
halo of mass M [= (47rr^/3)200pb, in terms of the virial 
radius r v ] found in CDM simulations O] 



p(r, M) 
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where /(c) = c 3 /[ln(l + c) — c/(l + c)] and the concen- 
tration parameter c depends weakly on mass. This profile 
implies an r _1 cusp for r < r v /c which will be altered 
by the presence of the Jeans scale. Solving for the Jeans 
scale in the halo rjh as a function of its mass using the 
enclosed mean density yields 

rj h ~ 3.4(c 10 //io) 1/3 m2 2 2/3 M 10 1/9 (n m / l 2 )- 2 / 9 kpc, (6) 

where we have scaled the mass dependent factors to the 
regime of interest c\q = c/10, /io = /(c)//(10), and 
Mio = M/10 10 Mq. For estimation purposes, we have 
assumed rjh <C r v /c which is technically violated for 
Mio ^, 1 and m = 10~ 22 eV, but with only a mild ef- 
fect. In the smallest halos, the Jeans scale is above the 
turnover radius r v /c, and there is no region where the 
density scales as r -1 . The maximum circular velocity 
will then be lower than that implied by eqn. (^). More 
massive halos will have their cuspy r~ l behavior extend 
from r — r v /c down to rjh- 

These simple scalings show that the wave nature of 
dark matter can prevent the formation of the kpc scale 
cusps and substructure in dark matter halos if m ~ 
10 _22 eV. However, alone they do not determine what 
does form instead. To answer this question, cosmologi- 
cal simulations will be required and this lies beyond the 
scope of the present work. Instead we provide here the 
tools necessary to perform such a study, a discussion of 
possible astrophysical implications and illustrative one 
dimensional simulations comparing FCDM and CDM. 

Linear Perturbations. — The evolution of fluctuations in 
the linear regime provides the initial conditions for cos- 
mological simulations and also directly affects the abun- 
dance of dark matter halos. Because the initial conditions 
are set while the fluctuations are outside the horizon, we 
must generalize the Newtonian treatment above to in- 
clude relativistic effects. 

Following the "generalized dark matter" (CDM) ap- 
proach of |L2[ ] , we remap the equations of motion for the 
the scalar field in equation (fil) onto the continuity and 



Euler equations of a relativistic imperfect fluid. First, 
note that in the Newtonian approximation, the current 
density j cx -0*V-0 — ip^tp* plays the role of momentum 
density so that "probability conservation" becomes the 
continuity equation. The dynamical aspect of equation 
(0) then becomes the Euler equation for a fluid with an 
effective sound speed c 2 ff = k 2 /4a 2 m 2 , where k is the 
comoving wavenumber. 

This Newtonian relation breaks down below the Comp- 
ton scale which for any mode will occur when a < k/2m. 
In this regime, the scalar field is slowly-rolling in its po- 
tential rather than oscillating and it behaves like a fluid 
with an effective sound speed c 2 ff = 1 jl2). For our pur- 
poses, it suffices to simply join these asymptotic solutions 
and treat the FCDM as GDM with 



c 2 -I 1 ' 



a < k/2m , 
a > k/2m , 
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with no anisotropic stresses in linear theory. We have 
verified that the details of this matching have a negligi- 
ble effect on the results. Since the underlying treatment 
is relativistic, this prescription yields a consistent, covari- 
ant treatment of the dark matter inside and outside the 
horizon. The linear theory equations including radiation 
and baryons are then solved in the usual way but with 
initial curvature perturbations in the radiation and no 
perturbations in the FCDM (TJ]. 

The qualitative features of the solutions are easily un- 
derstood. The comoving Jeans wavenumber scales with 
the expansion as kj cx ap 1 J i {a) or oc a 1 / 4 during matter 
domination and constant during radiation domination. 
Because the comoving Jeans scale is nearly constant, per- 
turbation growth above this scale generates a sharp break 
in the spectrum. More precisely, the critical scale is kj 

1 /2 — I 

at matter-radiation equality fcj cq = 9 vn^ 

' Mpc 1 . Nu- 
merically, we find that the linear density power spectrum 
of FCDM is suppressed relative to the CDM case by 



FCDM 



(k) = T 2 (fc)P C DM(fe) , T F (k) 
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where x = 1.61 m^^fc/fcjcq- The power drops by a fac- 
tor of 2 at 
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4.5m% 9 Mpc 
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The break in k is much sharper than those expected from 
inflation [|| or quartic self-interaction of a scalar field. In 
the latter, the Jeans scale is fixed in physical coordinates 
so that the suppression is spread over 3-4 orders of mag- 
nitude in scale Q]. 

Low Mass Halos. — In the CDM model, the abundance 
of low mass halos is too high when compared with the 
luminosity function of dwarf galaxies in the Local Group 
|l4],[l5|]. Based on analytic scalings, Kamionkowski & Lid- 
die T2|j argued that a sharp cutoff in the initial power 
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spectrum at k = 4.5ft- Mpc might solve this problem. 
Thus, the FCDM cutoff at k ~ 4.5 Mpc" 1 , produced 
if m ~ 10 _22 eV is chosen to remove kpc scale cusps, 
may solve the low mass halo abundance problem as well. 
Whether the required masses actually coincide in detail 
can only be addressed by simulations. 

Numerical simulations of CDM with a smooth cutoff 
in the initial power spectrum qualitatively confirm the 
analytic estimates but suggest that a somewhat larger 
scale may be necessary: half power at k = 2/iMpc" 1 
reduces the z = 3 abundance of 10 10 ft -1 M Q halos by a 
factor of ~ 5 and the abundance of 10 11 h~ x M Q halos 
by a factor of ~ 3 Note, however, that our model 

produces a much sharper cutoff in the power spectrum 
than in the model tested. Furthermore, astrophysical 
influences such as feedback or photoionization may have 
prevented dwarf galaxy halos from accumulating much 
gas or stars |l7j . 

First Objects and Reionization. — At very high redshift, 
much of the star formation in a CDM model is pre- 
dicted to occur in low-mass halos which are not present 
in the FCDM model. In a CDM model the first round 
of star formation is thought to occur in objects of mass 
~ 10 5 M Q (see Q and references therein) due to molec- 
ular hydrogen cooling. The consequent destruction of 
molecular hydrogen 19 2(}] implies that it is larger mass 
objects > 10 s M Q , where atomic cooling is possible, that 
are responsible for reionization. In our scenario, if the 
cutoff scale in eqn. (^) were set to reduce the abundance 
of M < 10 9 Mq halos, reionization could be delayed and 
the number of detectable galaxies prior to reionization 
reduced by a factor of ~ 5 ^jj . 

Live Halos. — Precisely what effect the Jeans (de Broglie) 
scale has on the structure and abundance of low mass 
halos is best answered through simulations. To provide 
some insight on these issues we conclude with simulations 
of the effects in one dimension. 

We solve the wave equation (3) in an interval < x < 
L with boundary ^(0) = ip(L) =0. At t = 0, the density 
perturbation is Sp = pg sm(Ttx/L) 3> Pb, with %j) real. 
We define the Jeans length rj by equation (4) with the 
density po] then the choice of rj/L specifies to. It is 
also convenient to define the dynamical time scale idyn = 
(47rGpo) -1 ^ 2 - A one dimensional CDM simulation with 
the same initial density and zero initial velocity was run 
for comparison. 

For rj ^> L, the field model does not form a gravi- 
tating halo. For rj ~ L a gravitationally bound halo is 
formed, but the cusp, which is clearly seen in the CDM 
simulation, is not observed (see Fig. [I]). Interference ef- 
fects cause continuous evolution on the dynamical time 
scale £dyn (cf. H). Note however that the gravitational 
acceleration is much smoother so that the trajectories of 
test particles (i.e. visible matter) will be less affected by 
these wiggles. 
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FIG. 1. One dimensional simulations (a) large Jeans scale 
rj/L = 0.44 (b) small Jeans scale rj/L = 0.044. Two snap- 
shots, i/tdyn = 99 (solid) and 100 (dotted) are shown for the 
density profile (units of 15po) and gravitational acceleration 
(units of 3Z//2tj yn , offset for clarity). 



For rj <C L, the density computed from the field fol- 
lows the CDM simulation when smoothed over many 
Jeans scales. This is to be expected because in this limit 
the Schrodinger equation can be solved in the geometrical 
optics approximation. Nonetheless, interference features 
localized in space (dx ~ rj) and time (5t < idyn) are 
quite strong, of order 1. Again these small-scale features 
make only a small contribution to the gravitational ac- 
celeration. 

Discussion. — We have shown that the wave properties of 
ultra light (m ~ 10 _22 eV) dark matter can suppress kpc 
scale cusps in dark matter halos and reduce the abun- 
dance of low mass halos. While such a mass may seem 
unnatural from a particle physics standpoint |p2| , even 
lighter scalar fields (to < 10~ 33 eV, where the wave na- 
ture is manifest across the whole particle horizon) have 
been proposed to provide a smooth energy component to 
explain observations of accelerated expansion p3] . 
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Three dimensional numerical simulations are required 
to determine whether our proposal works in detail. Our 
one dimensional simulations suggest that the small-scale 
cutoff appears at r ~ rj and that the density profile 
on these scales not only is not universal but also evolves 
continuously on the dynamical time scale (or faster) due 
to interference effects. The observable rotation curves 
are smoother than the density profile so that this pre- 
diction, while testable with high-resolution data, is not 
obviously in conflict with the data today. Likewise, the 
time-variation of the potential is smaller than that of 
the density but can in principle transfer energy from the 
FCDM to the baryons in the halo. This could puff up 
the baryons in dwarf galaxies while bringing the FCDM 
closer to a stationary ground state, but the precise evo- 
lution requires detailed calculations. 

Our Jeans scale is a weak function of density, rj oc 
p-l/i^ This has two testable consequences. The first is a 
sharp cutoff at k ~ 4.5 m 22 Mpc in the linear power. 
Quantities related to the abundance of low mass halos, 
e.g., dwarf galaxies in the Local Group, the first objects, 
faint galaxies at very high redshifts, and reionization can 
be seriously affected by the cutoff. Counterintuitively, 
quantities related to the non-linear power spectrum of 
the dark matter are only weakly affected due to the grav- 
itational regeneration of small-scale power []l6| . The sec- 
ond consequence is that choosing the mass to set the core 
radius in one class of dark matter dominated objects sets 
the core radius in another set, given the ratio of charac- 
teristic densities. This relation can in principle be tested 
by comparing the local dwarf spheroidals p4| to higher 
mass systems. 

While the detailed implications remain to be worked 
out, fuzzy cold dark matter provides an interesting means 
of suppressing the excess small-scale power that plagues 
the cold dark matter scenario. 
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